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Spin polarized states in dense neutron matter with BSk20 Skyrme force are considered in magnetic fields up to 10 G. 
It is shown that the appearance of the longitudinal instability in a strong magnetic field prevents the formation of a 
fully spin polarized state in neutron matter, and only the states with moderate spin polarization can be developed. 
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1. INTRODUCTION. BASIC EQUATIONS 

Magnetars are strongly magnetized neutron stars 
with emissions powered by the dissipation of mag- 
netic energy. The magnetic field strength at the sur- 
face of a magnetar is of about 10 14 -10 15 G. Such huge 
magnetic fields can be inferred from observations of 
magnetar periods and spin-down rates, or from hy- 
drogen spectral lines. In the interior of a magne- 
tar the magnetic field strength could reach values up 
to 10 20 G [1]. Then the issue of interest is the be- 
havior of neutron star matter, which further will be 
approximated by pure neutron matter, in a strong 
magnetic field |2J. In particular, a scenario is pos- 
sible in which a field-induced ferromagnetic phase 
transition occurs in the magnetar core. This idea 
was explored in the recent research [3], where it was 
shown that a fully spin polarized state in neutron 
matter could be formed in the magnetic field larger 
than 10 19 G. Note, however, that the breaking of the 
0(3) rotational symmetry in such ultrastrong mag- 
netic fields results in the anisotropy of the total pres- 
sure, having a smaller value along than perpendic- 
ular to the field direction [I] 5]. The possible out- 
come could be the gravitational collapse of a mag- 
netar along the magnetic field, if the magnetic field 
strength is large enough. Thus, exploring the possi- 
bility of a field-induced ferromagnetic phase transi- 
tion in neutron matter in a strong magnetic field, the 
effect of the pressure anisotropy has to be taken into 
account because this kind of instability could prevent 
the formation of a fully polarized state in neutron 
matter. In the present study, we determine thermo- 
dynamic quantities of strongly magnetized neutron 
matter taking into account this effect. 



Let us stop on the basic equations of the theory. 
The normal (nonsupcrfluid) states of neutron matter 
are described by the normal distribution function of 
neutrons f KlK2 = Tr 0O+a rei , where k = (p, a), p is 
momentum, a is the projection of spin on the third 
axis, and g is the density matrix of the system [SJ E] . 
The energy of the system is specified as a functional 
of the distribution function f,E = E(f), and deter- 
mines the single particle energy [7J [3] 

@jK2Kl 

The self-consistent matrix equation for determining 
the distribution function / follows from the minimum 
condition of the thermodynamic potential [7] and is 

/ = {exp(Y £ + Yt ■ Uncn + Y A ) + ly 1 (2) 
= {expCToO + IF 1 • 

Here the quantities e, Y, and I4 are matrices in the 
space of k variables, with (1^,4) k = Yi^5 KlK2 , 
Y = 1/T, Y l = -Hi/T and Y 4 = -fi /T being the 
Lagrange multipliers, [1q being the chemical poten- 
tial of neutrons, and T the temperature. In Eq. {2}, 
jtt„ = — 1.9130427(5)/iAr is the neutron magnetic mo- 
ment (fiN being the nuclear magneton), <jj are the 
Pauli matrices. 

Further it will be assumed that the third axis is 
directed along the external magnetic field H. Given 
the possibility for alignment of neutron spins along or 
opposite to the magnetic field H, the normal distri- 
bution function of neutrons and the matrix quantity 
£ (which we will also call a single particle energy) can 
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be expanded in the Pauli matrices <Xj in spin space 



/(P) = /o(p)ct 
£(p) = £o(p)oo 



/3(p)ct3, 
6(P)^3- 



(3) 
(4) 



The distribution functions /o, /a satisfy the nor- 
malization conditions 



v5>(p) = <?, (5) 
p 

X! /s(p) = £t - ej. = Ag. (6) 



Here g = g^ + is the total density of neutron mat- 
ter, g-f and g± are the neutron number densities with 
spin up and spin down, respectively. The quantity 
Ag may be regarded as the neutron spin order pa- 
rameter which determines the magnetization of the 
system M — /i n Ag. The magnetization may con- 
tribute to the internal magnetic field B = H + 47rM. 
However, we will assume, analogously to the previ- 
ous studies 2 , that, because of the tiny value of 
the neutron magnetic moment, the contribution of 
the magnetization to the inner magnetic field B re- 
mains small for all relevant densities and magnetic 
field strengths, and, hence, B w H. In order to get 
the self-consistent equations for the components of 
the single particle energy, one has to set the energy 
functional of the system. It represents the sum of the 
matter and field energy contributions 



E(f,H) = E m (f) + ^-V. 



(7) 



The matter energy is the sum of the kinetic and 
Fermi-liquid interaction energy terms [5J |S] 



E m (f) = E (f) + E int (f), 

£ (/) = 2^e (p)/o(p), 
p 

Eint(f) = X){eo(p)/o(p) + ea(p)/s(p)}, 



(8) 



where 



£o(p) = ^ E ^"(k)/o(q), k = 5—5, 
q 

£ 3 (p) = ^E c/ "( k )/3(q)- 



(9) 
(10) 



Here e (p) = is the free single particle spectrum, 
to is the bare mass of a neutron, Uq (k), f/f(k) are 
the normal Fermi liquid (FL) amplitudes, and £o,£3 
are the FL corrections to the free single particle spec- 
trum. Taking into account Eqs. Ql,© and ©, ex- 
pressions for the components of the single particle 
energy read 

£o(p) = £o(p) + £q(p) - Mo, &(p) = -pi n H + £3{p)- 

(11) 



In Eqs. the quantities £~o, £3 are the function- 
als of the distribution functions fa , which, using 
Eqs. and (J3j> , can be expressed, in turn, through 
the quantities £: 



fo = ^{n(u}+) + n(w_)}, 
h = 7}{n(u+) - n(u)-)}, 



(12) 
(13) 



where 



n(u±) = {exp(Y uj±) + 1} 1 , u± = £ ± f 3 . 

Thus, Eqs. (fTTj) - ([T3"|) form the self-consistency 
equations for the components of the single particle 
energy, which should be solved jointly with the nor- 
malization conditions ([5]), 

The pressures (longitudinal and transverse with 
respect to the direction of the magnetic field) in the 
system are related to the diagonal elements of the 
stress tensor whose explicit expression reads [5] 



O'ik 



f- 



8f 

do 



H,T 



HiBk 

47T 



Here 



f=f H - 



47T ' 



(14) 



(15) 



f H = jj(E- TS) - HM is the Helmholtz free energy 
density. For the isotropic medium, the stress ten- 
sor (|14p is symmetric. The transverse pt and longitu- 
dinal pi pressures are determined from the formulas 



Pt 



-011 



-022, Pi 



"0'33- 



At zero temperature, using Eqs. (JZJ, 
get the approximate expressions 



Pt = g\ 



Pi = Q\ 





\ H2 


dg / 






\ H2 


dg / 


>h &m Stt' 



(16) 
(17) 



where e m is the matter energy density, and we disre- 
garded the terms proportional to M. In ultrastrong 
magnetic fields, the quadratic on the magnetic field 
term (the Maxwell term) will be dominating, leading 
to increasing the transverse pressure and to decreas- 
ing the longitudinal pressure. Hence, at some critical 
magnetic field, the longitudinal pressure vanishes, re- 
sulting in the longitudinal instability of neutron mat- 
ter. The question then is: What is the magnitude of 
the critical field and the corresponding maximum de- 
gree of spin polarization in neutron matter? 



2. EOS OF DENSE NEUTRON MATTER 
IN A STRONG MAGNETIC FIELD 

In numerical calculations, we utilize the BSk20 
Skyrme force |10) constrained such as to avoid the 
spontaneous spin instability of neutron matter at 
densities beyond the nuclear saturation density and 
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to reproduce a microscopic EoS of nonpolarized neu- 
tron matter. Expressions for the normal FL ampli- 
tudes in Eqs. (l9|). (fT0|) in terms of the parameters of 
the Skyrme interaction are given in Ref. 111] . Now 
we present the results of the numerical solution of 
the self-consistency equations. Fig. 1 shows the spin 
polarization parameter II = ^ of neutron matter 
as a function of the magnetic field H at two differ- 
ent values of the neutron matter density, g = 3 go 
and g = Ago, which can be relevant for the magnetar 
core. It is seen that the impact of the magnetic field 
remains small up to the field strength 10 17 G. The 
larger the density is, the smaller the effect produced 
by the magnetic field on spin polarization of neutron 
matter. 
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Fig.l. Neutron spin polarization parameter as a 
function of the magnetic field H for the Skyrme 
force BSk20 at zero temperature and fixed values 
of the density, g = 3g<j ond g = 4go- The vertical 
arrows indicate the maximum magnitude of spin 
polarization attainable at the given density, see 
further details in the text. 



At the magnetic field H = 10 18 G, usually con- 
sidered as the maximum magnetic field strength in 
the core of a magnetar (according to a scalar virial 
theorem, see Ref. [T] and references therein) , the mag- 
nitude of the spin polarization parameter doesn't ex- 
ceed 33% at g = 3go and 18% at q = 4^ - However, 
the situation changes if the larger magnetic fields are 
allowable: With further increasing the magnetic field 
strength, the magnitude of the spin polarization pa- 
rameter increases till it reaches the limiting value 
II = — 1, corresponding to a fully spin polarized state. 
For example, this happens at H m 1.25 • 10 19 G for 
g = 3go and at H w 1.98 • 10 19 G for g = 4g , i.e., 
certainly, for magnetic fields larger than 10 19 G. Nev- 
ertheless, we should check whether the formation of a 
fully spin polarized state in a strong magnetic field is 
actually possible by calculating the anisotropic pres- 
sure in dense neutron matter. The meaning of the 
vertical arrows in Fig. 1 is explained later in the text. 
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Fig. 2. (a) Pressures, longitudinal ( descending 
branches) and transverse (ascending branches), as 
functions of the magnetic field H for the Skyrme 
force BSk20 at zero temperature and fixed values 
of the density, g = 3g$ and g = 4go- (b) Same as 
in the top panel but for the normalized difference 
between the transverse and longitudinal pressures. 

Fig. 2a shows the pressures (longitudinal and 
transverse) in neutron matter as functions of the 
magnetic field H at the same densities, g = 3go 
and g = Ago. First, it is clearly seen that up to 
some threshold magnetic field the difference between 
transverse and longitudinal pressures is unessential 
that corresponds to the isotropic regime. Beyond 
this threshold magnetic field strength, the anisotropic 
regime holds for which the transverse pressure in- 
creases with H while the longitudinal pressure de- 
creases. The longitudinal pressure vanishes at some 
critical magnetic field H c marking the onset of the 
longitudinal collapse of a neutron star. For example, 
H c w 1.56 • 10 18 G at g = 3g and H c w 2.42 ■ 10 18 G 
at g = Ago. In all cases under consideration, this 
critical value doesn't exceed 10 19 G. 

The magnitude of the spin polarization parameter 
II cannot also exceed some limiting value correspond- 
ing to the critical field H c . These maximum values 
of the II's magnitude are shown in Fig. 1 by the ver- 
tical arrows. In particular, II C « —0.46 at g = 3go 
and II C w —0.38 at g = 4go- As can be inferred from 
these values, the appearance of the negative longitu- 
dinal pressure in an ultrastrong magnetic field pre- 
vents the formation of a fully spin polarized state in 
the core of a magnetar. Therefore, only the onset 
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of a field-induced ferromagnetic phase transition, or 
its near vicinity, can be catched under increasing the 
magnetic field strength in dense neutron matter. A 
complete spin polarization in the magnetar core is 
not allowed by the appearance of the negative pres- 
sure along the direction of the magnetic field, con- 
trary to the conclusion of Rcf . [3] where the pressure 
anisotropy in a strong magnetic field was disregarded. 

Fig. 2b shows the difference between the trans- 
verse and longitudinal pressures normalized to the 
value of the pressure po in the isotropic regime (which 
corresponds to the weak field limit with pi = p t — po) 
being 8 = Pt ~ Pl ■ Applying for the transition from the 
isotropic regime to the anisotropic one the criterion 
8 ~ 1, the transition occurs at the threshold field 
H th « 1.15-10 18 G for g = 3g and H th w 1.83-10 18 G 
for g = 4qq. In all cases under consideration, the 
threshold field H t h is larger than 10 18 G, and, hence, 
the isotropic regime holds for the fields up to 10 18 G. 
The vertical arrows in Fig. 2b indicate the points cor- 
responding to the onset of the longitudinal instability 
in neutron matter. The maximum allowable normal- 
ized splitting of the pressures corresponding to the 
critical field H r is 8 ~ 2. 
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Fig. 3. Same as in Fig. 2 but for: (a) the Helmholtz 
free energy density of the system; (b ) the ratio of the 
magnetic field energy density to the Helmholtz free 
energy density of the system. 

Fig. 3a shows the Helmholtz free energy density of 
the system as a function of the magnetic field H. It is 
seen that the magnetic fields up to H ~ 10 18 G have 
practically small effect on the Helmholtz free energy 



density f#, but beyond this field strength the contri- 
bution of the magnetic field energy to the free energy 
fjy rapidly increases with H. However, this increase 
is limited by the values of the critical magnetic field 
corresponding to the onset of the longitudinal insta- 
bility in neutron matter. The respective points on 
the curves are indicated by the vertical arrows. 
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Fig. 4- The Helmholtz free energy density of the 
system as a function of: (a) the transverse pressure 
Pt, (b) the longitudinal pressure pi for the Skyrme 
force BSk20 at zero temperature and fixed values of 
the density, g — 3 go and g = 4go- 

Fig. 3b shows the ratio of the magnetic field en- 
crgy density e/ = to the Helmholtz free energy 
density at the same assumptions as in Fig. 2. The in- 
tersection points of the respective curves in this panel 
with the line et/fji = 0.5 correspond to the magnetic 
fields at which the matter and field contributions to 
the Helmholtz free energy density are equal. This 
happens at H w 1.18 • 10 18 G for g = 3go, and at 
H w 1.81 • 10 18 G for g = Ag Q . These values are 
quite close to the respective values of the threshold 
field H t h, and, hence, the transition to the anisotropic 
regime occurs at the magnetic field strength at which 
the field and matter contributions to the Helmholtz 
free energy density become equally important. It is 
also seen from Fig. 3b that in all cases when the lon- 
gitudinal instability occurs in the magnetic field H c 
the contribution of the magnetic field energy density 
to the Helmholtz free energy density of the system 
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dominates over the matter contribution. 

Because of the pressure anisotropy, the EoS of 
neutron matter in a strong magnetic field is also 
anisotropic. Fig. 4 shows the dependence of the 
Helmholtz free energy density ijj of the system on 
the transverse pressure (top panel) and on the longi- 
tudinal pressure (bottom panel) at the same densities 
considered above. Since in an ultrastrong magnetic 
field the dominant Maxwell term enters the pressure 
p t and free energy density fjj with positive sign and 
the pressure pi with negative sign, the free energy 
density f# is the increasing function of p t and de- 
creasing function of pi. In the bottom panel, the 
physical region corresponds to the positive values of 
the longitudinal pressure. 

The obtained results can be of importance in the 
structure studies of magnetars. It would be also of 
interest to extend this research to finite tempera- 
tures relevant for proto-neutron stars which can lead 
to a number of interesting effects, such as, e.g., an 
unusual behavior of the entropy of a spin polarized 
state QUE?]. 

J.Y. was supported by grant 2010-0011378 from 
Basic Science Research Program through NRF of Ko- 
rea funded by MEST and by grant R32-10130 from 
WCU project of MEST and NRF. 

References 

1. E. J. Ferrer, V. de la Incera, J. P. Keith, I. Por- 
tillo, and P. L. Springsteen. Equation of state of 
a dense and magnetized fermion system / / Phys. 
Rev. C 2010, 82, 065802, p. 15. 

2. A. A. Isayev, and J. Yang. Spin-polarized states 
in neutron matter in a strong magnetic field / / 
Phys. Rev. C 2009, 80, 065801, p. 7. 

3. G. H. Bordbar, Z. Rezaei, and A. Montakhab. 
Investigation of the field-induced ferromagnetic 
phase transition in spin-polarized neutron mat- 
ter: A lowest order constrained variational ap- 
proach // Phys. Rev. C 2011, 83, 044310, p. 7. 



4. V. R. Khalilov. Macroscopic effects in cold mag- 
netized nucleons and electrons with anomalous 
magnetic moments / / Phys. Rev. D 2002, 65, 
056001, p. 6. 

5. A. A. Isayev, and J. Yang. Spin polarized states 
in strongly asymmetric nuclear matter / / Phys. 
Rev. C 2004, 69, 025801, p. 8. 

6. A. A. Isayev. Spin ordered phase transitions in 
isospin asymmetric nuclear matter / / Phys. Rev. 
C 2006, 74, 057301, p. 4. 

7. A. I. Akhiezer, A. A. Isayev, S. V. Peletminsky, 
A. P. Rekalo, and A. A. Yatsenko. On a theory 
of superfluidity of nuclear matter based on the 
Fermi-liquid approach // JETP 1997, 85, 1-12. 

8. A. I. Akhiezer, A. A. Isayev, S. V. Peletmin- 
sky, and A. A. Yatsenko. Multi-gap superfluid- 
ity in nuclear matter // Phys. Lett. B 1999, 451, 
430-436. 

9. L. D. Landau, E. M. Lifshitz, and L. P. 
Pitaevskii. Electrodynamics of Continuous Me- 
dia, New York, Pergamon, 1984, 2nd ed. 

10. S. Goriely, N. Chamel, and J. M. Pearson. 
Further explorations of Skyrme-Hartree-Fock- 
Bogoliubov mass formulas. XII. Stiffness and sta- 
bility of neutron-star matter / / Phys. Rev. C 
2010, 82, 035804, p. 18. 

11. A. A. Isayev, and J. Yang. Phase transition to the 
state with nonzero average helicity in dense neu- 
tron matter // JETP Lett. 2010, 92, p. 783-787. 

12. A. A. Isayev. Finite temperature effects in anti- 
ferromagnetism of nuclear matter / / Phys. Rev. 
C 2005, 72, 014313, p. 4. 

13. A. A. Isayev. Unusual temperature behavior of 
the entropy of the antifcrromagnetic spin state 
in nuclear matter with an effective finite range 
interaction // Phys. Rev. C 2007, 76, 047305, 
p. 4. 



5 



